ob ■ The Darboux process and a 

§ ; noncommutative bispectral problem 

(N 

F. Alberto Griinbaum 

I Abstract. The Darboux process, also known by many other names, 

■ played a very important role in some extremely enjoyable joint 
I work that Hans and I did 25 years ago. I revisit a version of this 

problem in a case when scalars are replaced by matrices, i.e., ele- 
p I I ments of a non-commutative ring. Many of the issues studied here 

■ can be pushed to the case of a ring with identity, but my emphasis 
(-H I is on very concrete examples involving 2x2 matrices. 

» 

a ■ 

The material in this paper was presented at a conference in honor of 
Hans Duistermaat, Aug 20-24 2007. 

t:J- ' 1. The bispectral problem 

m ■ 

OO ■ About 30 years ago I was extremely lucky. 1 gave a talk in Berkeley 

where I mentioned the following problem: 

Find all nontrivial instances where a function ip{x, k) satisfies 



(N 

od 
O . 

O ■ ( ^' j~ ) ^i-"^-! ^) = (— -D^ + V{x))(f{x, k) = k(p{x, k) 



>< 



dx 



as well as 



^. : B (^k, ^{x, k) ^ [y,^ h{k) j ^(x, k) = e{x)^{x, k). 

All the functions ^(x), 9(x) are, in principle, arbitrary except 
for smoothness assumptions. Notice that here M is arbitrary (finite). 

1 was fortunate that Hans was in the audience, and about a week 
later he came up with a tool to attack this problem. After a few months 
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of intense work, mainly by slow mail, we found ourselves with a rather 
nice picture. 

The complete answer to the problem is given as follows: 

Theorem 1 |DG] . If M = 2, then V{x) is (except for translation) 
either c/x^ or ax, i.e., we have a Bessel or an Airy case. If M > 2, 
there are two families of solutions. 

a) L is obtained from Lq = —D"^ by a finite number of Darboux 
transformations {L = AA* — > L = A* A). In this case V is 
a rational solution of the Korteweg-deVries equation and all 
rational solutions of KdV decaying at infinity show up in this 
fashion. 

b) L is obtained from Lq = —D^ + after a finite number of 
rational Darboux transformations. 

It was later observed in [MZ j that in the second case we are dealing 
with rational solutions of the Virasoro or master symmetries of KdV. 

In the first case the space of common solutions has dimension one, 
in the second it has dimension two. One refers to these as the rank one 
and rank two situations. In jPG ] one finds several other equivalent 
descriptions of the solution such as those in terms of the monodromy 
group of the equation. 

Observe that the "trivial cases" when M = 2 are self-dual in the 
sense that since the eigenfunctions f{x, k) are functions either of the 
product xk or of the sum x + k, one gets B by replacing k for x in L. 

My reasons for asking this question can be traced back to an effort 
to understand some work on "time-and-band-limiting" that had led 
me to isolate certain properties of well known special functions. For 
an example relating to orthogonal polynomials see |Glj . For a more 
updated versions of this connection between the bispectral problem and 
the issue of time-and-band-limiting see, |G2L IG31 IGY2j . 

The work with Hans gave rise to a large number of papers by other 
people, some of which can be found in the list at the end of this paper. 
This listing is far from complete and I apologize for the omissions. 

It may be appropriate to observe that what we are calling the Dar- 
boux process has been reinvented many times, including in the work 
of some rather well known people, see for instance [Sc|, lIHj . Reference 
jYZ] talks about the Geronimus transformation, from 1940, and its in- 
verse the Christoffel transform. It is clear that the first one (as noticed 
in |YZ] ) has a lot in common with what we are calling the Darboux 
transformation. See also |SVZL [Z] . 
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2. The Bochner Krall problem 

In a series of papers with Luc Haine, [GHIL IGH2L IGH3] and then 
with Luc and Emil Horozov, see [GHHIL IGHH2j . we noticed that a 
large class of polynomials 

Pn{k) 

that satisfy three term recursion relations in the variable n, as well as 
differential equations in the variable k can be obtained by an applica- 
tion of a similar Darboux transformation starting from the so called 
classical orthogonal polynomials of Jacobi, Laguerre and Hermite. In 
this case one goes from a tridiagonal matrix 

Lo = AB 

to a new tridiagonal matrix L by factorizing the first one (or a function 
of it) as a product of two bidiagonal matrices. As indicated in |GHlj 
some form of this method is given in [MSlj . 

In this case one runs into the Toda fiows and its master symmetries. 
Further work on these lines can be found in [GYlj and |GH4L IGH5j 
and for a very nice survey of all this material see |H2j . 

The origins of this line of work is contained in papers such as |BoL 

3. A matrix valued version of the Darboux process for a 

difference operator 

Consider the block tridiagonal matrix Lq 

/Bo I 

where all the matrices Ai, Bi are of size N x N and / denotes the N xN 
identity matrix. 

If we try to factorize this in the form 

Lq = a[3 

where 

a 



and 



ai / 
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with all the matrices Oj, /3j are of size NxN, and then define the matrix 

L = (3a 

we have that 

[Bo I \ 

L= Ai B, I 

\ ■••/ 

where all the matrices Ai,Bi are of size NxN, and the following 
relations hold 

which then gives 

Bn = Bn — Pn+l + = -^n-l — Q^n-l + «n 

and 

An = /3nAn-l/3j^_i = AnQn-l- 

These expressions are valid for n = 2, 3, . . . in the case of An and for 
n = 1,2, . . . in the case of Bn- 

Above we take Po = 0, so that Bq = = Bq — Pi. We also need 
to take Ai = {Bq — ao)aQ. 

A moment's thought gives that once Lq is given, the only free pa- 
rameter is the matrix a^. 

Just as in [GHlj . and in spite of the fact that one is dealing here 
with a semi-infinite block tridiagonal matrix, it is possible to see the 
connection between this construction and that in [MSlj . One puts 

Pn = -4>n(f)n-l 

and then notices that this amounts to choosing (f) in the null-space of 
L. Since L is not doubly- infinite we seem to have lost some freedom 
in picking this subspace, but this can be remedied as in [GHl] by 
considering L as a limit of an appropriate doubly-infinite matrix with 
a rich null-space. 

4. Fancier versions of the Darboux process 

It is well known that it is useful to extend the original method of 
Darboux consisting of going from 

Lo = AB 

to 

L = BA 

in an appropriate way. Notice that in the standard case we have 

BLo = LB. 
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The new idea is to allow for an arbitrary banded matrix ( or a 
differential operator) U and to declare L a Darboux transform of Lq as 
long as we have 



In several of the uses of the original Darboux's method one needs 
to apply it repeteadly and this fancier version of the method takes care 
of that. 

One should also keep in mind the results in [GHHIL [GHH2j IH2j 

where the usual factorization followed by a reversal of the factors is 
applied not directly to L but to a constant coefficient polynomial in L. 

I thank Jose Liberati for pointing out to me that at the very end 
in [GGRW] one finds an application of the theory of quasidetermi- 
nants (a notion that goes back to Cayley) to obtain expressions for 
matrix valued orthogonal polynomials in terms of their matrix valued 
moments. Many of these results, as well as others, have been derived 
independently by making good use of the notion of Schur complements 
in L. Miranian's Berkeley thesis. The main new results are contained 
in [Mj . In the next section we give a very brief look into the the- 
ory of matrix valued orthogonal polynomials and a short guide to the 
literature that is relevant to us. 

One should remark that this same theory of quasideterminants has 
been studied in connection with a certain Darboux process for a matrix 
Schroedinger equation in |GV] . In this case one needs to consider this 
fancier version. For a nice use of quasideterminants in our context see 



The matrix version of the Darboux process for a difference opera- 
tor discussed in the previous section could be extended in this fancier 
fashion too. 

5. Matrix valued orthogonal polynomials 

Given a self-adjoint positive definite matrix valued weight function 
W{t), Krein, see |Kll IK2] . considers the skew symmetric bilinear form 
defined for any pair of matrix valued functions P{t) and Q{t) by the 
numerical matrix 



where Q*{t) denotes the conjugate transpose of Q{t). 

Proceeding as in the case of a scalar valued inner product Krein 
proves that there exists a sequence {Pn)n of matrix polynomials, or- 
thogonal with respect to W, with P„ of degree n and monic. 



LU. 



m 
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Krein goes on to prove that any sequence of monic orthogonal ma- 
trix valued polynomials {Pn)n satisfies a three term recurrence relation 

(1) ^nPn-l(t) + 5„P„(t) + P„+i(t) = tP„(t), 

where P_i is the zero matrix and Pq is the identity matrix. These 
coefficient matrices enjoy certain properties: in particular the An are 
nonsingular. 

The equation above can be rewritten as 

with a matrix ^ such as the one that has appeared in previous sections. 

To place ourselves in the context of the bispectral problem we 
consider matrix valued polynomials (P„)„ satisfying not only the equa- 
tion above but also "right hand side" differential equations of the form: 

(2) P„D = A„P„ for all n>0 

with A„ a matrix valued eigenvalue and D a differential operator of 
order s with matrix coefficients given by 

i=0 

which acts on P{x) by means of 

s 

PD = J29"iP){x)Fi{x). 

1=0 

This problem in the matrix case has been studied in [D], IDGIL 
[GToL [GIL IGPTIL IGPT2L IGPT3L IGT] and in a few other places. 

One can see that the differential operators that correspond to a 
fixed family of polynomials form an associative algebra which in gen- 
eral is non- commutative, see |DG1L IGTj . The problem of exhibiting 
elements of this algebra that have a minimal order will occupy us in a 
few examples in the next two sections. 

6. A few examples 

Here we consider in detail few examples of the matrix version of the 
basic Darboux process described above. 

For A > 3/2 consider the monic matrix valued polynomials which 
are orthogonal with respect to the weight matrix 

W{x) = ((2 - x)xr--'/' ""-^y xe[0,2]. 
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Let 

'Bo I 



be the corresponding block tridiagonal matrix with 

2(n + A)(n + A- 1) 
_ n(n + 2A - 2) 
" " 4(n + A-l)2^- 

Here S = l] and / ' 



1 J \0 1^ 

These polynomials can be seen to be joint eigenfunctions of a first 
order differential operator, an observation that was made for a special 
value of A in [CGT1^G2] . 

If ao is an arbitrary matrix we can consider the monic polynomials 
that result from one application of the Darboux process to the block 
tridiagonal matrix Lq with free parameter ao- 

We can see that for an invertible symmetric the new orthogo- 
nality weight is given by 

1 x-1 
X - 1 1 



W{x) = (2 - x)^-3/2^A-5/2 

-22^5e(^,^) //2A-2 -1 



4(2A-3) VV -1 2A-2 



(2A - 3)ao ^ 6oix 



Here Be stands for the usual Beta function. 

Recall that the Darboux process played an important role in getting 
new bispectral situations in |DG] . 

We show below some examples that illustrate that for appropriate 
values of A the new polynomials are joint eigenfunctions of some higher 
order differential operators, i.e., we get new bispectral situations. This 
appears to have little to do with ao being symmetric. 

Example 1. 

A = 5/2 ao = (3 fj. 
Here we find one differential operator D satisfying 



PriD \fiPn 
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and 



A. 



D = J2 Br 



r=0 



S3 = A(t - 2)t(3t - 2) 

24 / 1 _i 

= yi(7t - 9) (^_\ / 

8 / 5i + 6 -8i 
^ ~ 5 V-2(5i + 3) m 

/ 8x11 _c 



(n+2) (5n3 +20n^+3n-44) 5n^ +30n3+43n^ - 14w+40 ' 

5 5 
5n^+30n3+43n^+2n+32 n(5n3+30n^+43n+26 



There are no operators of lower order in the algebra. 



Example 2. 



A = 7/2, ao = (^5 



Here the coresponding operator is given by 



6 
r=0 
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with 



15 \-l 1 
2(i - 2)2^2(5^ _ 4) / 1 _i 



{t-2)H'fl -1 

^ 5 
^4 = 4(i - 2)t{5t'^ -8t + 2) 

Bs = 16t{5t^ - m + 6) 

B2 



1 1 

1 -1 



-1 1 

1 -1' 
-1 1 



Bo 



and we have 



16t(131t- 148) / 1 -1 
39 V-1 1 

32 /-4(2t - 7) 3{t - 6) 
19 V (9* - 28) -2{2t - 9) 

J_ /-13 X 32 -3 X 64 
19 i 11 X 32 



with 



+ 



19n^ + 285n^ + 1615n^ + 42 7571^ + 2446/1^ 
285 

1 /-12480n - 6240 10080n - 2880 
285 V 12960n + 5280 -10560n 



1 



1 -1 
-1 1 



Example 3. 

A = 9/2 ao - ( Q Q I • 

In this case there is one operator of order eight whose corresponding 
A„ is given below 

(n - 3)(n + 6)(n + 7)(n + 8)an -{n - 2)n{n + 7)(n + 8)/?, ' 
-(n + l){n + 6)7„ (n - l)n{n + l)(n + 10)();, 



An = 

with 



tt„ = + lOn'"^ + 59n^ + 170n + 840 

/5„ = n^ + Un^ + 95^2 + 322n + 1080 

7„ = + 21n^ + 169n^ + 651n-'^ + 1198n^ + 840n - 20160 

5„ = + 18^3 + 14371^ + 558n +1512. 
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Once again, this corresponds to the lowest order possible differential 
operator in the corresponding algebra. 



7. A few Jacobi type examples 

A different avenue for exploring the similarities as well as the dif- 
ferences between the use of the Darboux process in the scalar and the 
matrix valued case is given by the examples in this section. 

First recall that in the scalar case it follows from results in [KoKolT, 
[Z], IH21 IGYlj that the polynomials orthogonal to the weight con- 
sisting of a Jacobi density plus two possible delta masses of nonnegative 
strengths W, V at the ends of the interval, i.e., 

= (1 - + xf + W6i{x) + V6.i{x) 

satisfy differential equations when a and f3 satisfy certain natural in- 
tegrality conditions. The simplest example is given by the so called 
Koorwinder polynomials, corresponding to a = f3 = 0. If the weight 
at 1 is the only one that is present then the order is 2a + 4. If both 
delta weights are thrown in, then the order is 2q; + 2/3 + 6. The results 
can be obtained by an application of the Darboux process as shown in 
jH2l [GYT] . 

We consider now a small colection of situations analogous to the 
ones above. 

The weight matrices will, as before, consist of a matrix weight den- 
sity plus a pair of deltas at the end points weighted by certain matrices 
W, V, i.e., we have 

W{x) = (1 - a;)"(l + xf Q + W5i{x) + V5.i{x). 

For the first batch of examples we will assume that a, (3 are both 0. If 
V and W coincide with the matrix ( ^ ] then we find two linearly 



^0 0^ 

independent operators of order 5 and one of order 6 as well as other 
operators of higher order. There are no other operators of lower order. 

If V is the matrix and W is the matrix or the ma- 

trix then we find two linearly independent operators of order 6 

as well as other operators of higher order. There are no other operators 
of lower order. 
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^2 ) S'lid W is the matrix 
^2 1 then we have the same situation as in the last example. 
In general if V and W are of the form ( ^ ^ j and ( ^ ] then 



b cj \e f 

the lowest order operator in the algebra is just one operator of order 8. 

We come now to a different sort of examples. 

Assume that a and P (> —1) are arbitrary, but insist in picking W 
and V to be arbitrary (and non-necessarily equal) nonnegative multi- 

p>esoftKen,atn.(; ; 

In this case there is a very nice second order differential operator 
in the algebra which is independent of the choice of the scalar factors 
that appear in front of the matrix above to give W and V. There is no 
lower order operator in the algebra. When the deltas are both missing 
then the algebra contains an operator of order 1. 

The right handed differential operator alluded to above is a scalar 
operator of the usual Jacobi type, with coefficients (1 — x"^) and {a + 

P — 1) — x{a + P + 3) multiplied on the right by the matrix ^ . 

The eigenvalue is —n{n + a + (3 + 2) multiplied by this same matrix. 

8. An explicit differential operator 

The paper [DG] contains a proof that in the continuous-continuous 
case, when all operators in question are differential operators, the so 
called ad-conditions 

adL'"+i(e) = 

are necessary and sufficient to have what has eventually been called 
a bispectral situation, i.e., a solution of the original problem tackled 
with Hans. 

This condition gives a set of nonlinear equations that need to be 
solved in the unknowns L, 0. 

It is important to see that this condition can be easily adapted 
to other situations, including the present noncommutative one. This 
approach was taken up in |GIj and in [GTj . In the second of these 
papers the "ad-conditions" are shown to be equivalent, once again, to 
bispectrality. 

In general finding the differential operators of lowest possible order 
that appear in a bispectral situation is not easy. By repeated appli- 
cations of the Darboux process one may end getting elements of the 
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corresponding algebra that are not necessarily of the lowest possible 
order. This issue has surfaced in several different papers, starting with 
|DG] and a nice account is given in [H2]. 

In |GTj one finds an explicit construction of a differential operator 
that results form the conditions 

adL'"+i(A) = 0. 

The operator D is given by 



D = ^d'{P)'- 



r=0 

with matrix coefficients Sk = Sk{x) given by 

Sk = {iL-xir-''AP)o. 
In particular, we display some of the coefficients. 
So = ((L - x/)™AP)o, 

and, at the other end, 

= (AP)o = Ao, 
and this operator satisfies the desired condition 

PnD = KPn n>0. 

The subindex above refers to the first entry of the corresponding 
"vector" with matrix valued entries. 

9. Toda flows with matrix valued time 

As seen in |GH2] . and certainly in other places too, repeated ap- 
plications of the scalar Darboux process introduces "times and flows" 
that are related to the Toda flows. Since these times appear as the free 
parameters in each application of the process, it is only natural to raise 
the issue of "matrix valued times" . 

10. Electrostatics: Heine, Stieltjes, Darboux 

In a remarkable paper that follows on earlier work of Heine, Steiltjes 
came up with a nice electrostatic interpretation for the zeros of the 
Jacobi polynomials. Later work of Stieltjes as well as of I. Schur and 
D. Hilbert showed similar interpretations in the case of the Laguerre 
and Hermite polynomials. 

In |G7L IG8j I raise the possibility of some relation between the Dar- 
boux process where the orthogonality functional gets more and more 
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complicated with every application of the process and the correspond- 
ing electrostatic interpretation of the families of polynomials that ap- 
pear along the way. 

It would be interesting to see what if anything of this picture can 
be developed in the matrix valued case. 

11. Markov chains 

In |G4L IG51 IG6|. IGdl] one finds examples of interesting quasi- 
birth-and-death processes that can be studied by exploting their con- 
nection with certain specific examples of matrix valued orthogonal 
polynomials. In particular in [G5L IG6] one finds examples where the 
recurrence of the process is related to the presence of a matrix val- 
ued delta weight at 1. Since the appearance of these delta weights is 
one of the main characteristic of an application of the Darboux pro- 
cess one may wonder about a probabilistic interpretation of the relation 
that may exist between two Markov chains whose transition probability 
matrices are related by a Darboux transformation. 

12. Things that appear before their time 

One of the most surprising phenomena uncovered in |DG] has to do 
with what was called "the cusps", namely degenerate situations that 
corespond to degeneracies of "higher order operators" yielding "lower 
order" ones. To put this in the context of scalar valued orthogonal 
polynomials consider the so called Koornwinder polynomials which are 
orthogonal to Lebesgue measure in [—1, 1] plus a pair of delta masses 
at the end points of the interval. In this case one knows that the 
corresponding orthogonal polynomials are the common eigenfunctions 
of a sixth order differential operator. 

In the special case when the strength of the two delta masses agrees 
one gets an operator of order four, and one can say that in a search 
according to the order of these operators this example, just as "the 
cusps" in |DG] . appears before its time. 

We made a tentative exploration of the situation in the matrix 
valued case and examples of this phenomenon are seen in section 7. 

13. The multivariable case 

In this section we mention that in [G9j one finds a specific random 
walk introduced by Hoare and Rahman, see [HRj which we show leads 
to a bispectral situation in terms of polynomials of two variables. I 
also want to mention that in the multivariable case one finds a version 
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of the Darboux process to obtain interesting deformations of the two 
dimensional Chebyshev measure, see |GIlj . 

14. Conclusion 

It is clear that very little of the development that I have tried to 
summarize here could have happened were it not for my good fortune 
of teaming up with Hans at the beginning of this journey. 
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